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ATMOSPHERIC TURBULENCE EFFECTS ON 
AIRCRAFT NOISE PROPAGATION 


Robert L. Chapkis 

INTRODUCTION 

Of all the characteristics of the atmosphere that affect sound propaga- 
tion, turbulence is the least understood . Attenuation of sound by viscous 
and heat conduction effects is now well understood and can be calculated 
(refs. 1, 2, 3). Refraction of sound rays caused by wind and temperature 
gradients can also be determined if profiles of the mean temperature and wind 
are known or measured (ref. 4). If the acoustical Impedance of the ground 
surface beneath the atmosphere is known or can be measured, then the sound 
attenuation caused by the ground cau be approximately calculated (ref. 5). 

Accounting for atmospheric turbulence effects, however, is much more 
difricult because there is still uncertainty as to the physical mechanisms 
by which turbulence affects sound propagation. Several theoretical studies 
have been based on a physical model which assures that the important aspect 
of turbulence is the scattering of an initially coherent sound wave by 
thermal and momentum fluctuations in the atmosphere (refs. 6 through 8). 

The loss in intensity of the coherent wave manifests itself as an apparent 
attenuation caused by turbulence. 

Other studies of turbulence effects on sound propagation assumed that 
the important effect is the broadening of a beam of sound caused by phase 
fluctuations induced by the turbulence (ref. 9). The increase in cross- 
sectional area of the beam causes a corresponding decrease in sound intensity 
because the sound-energy flux is spread out over a larger area. 

Finally, there have been other studies of sound attenuation due to 
turbulence which show that for certain conditions, such as for propagation 
of sound through a highly turbulent jet, actual turbulent absorption might 
be important (refs. 10 and 11). The absorption comes about through an 
interaction of a sound wave with a turbulent flow field involving a net loss 



of energy front Che sound wave. A process of cascading of energy from large 
to small turbulent eddies occurs whereby acoustic energy is ultimately 
dissipated through viscous forces acting on the smallest eddies. 

In addition to the various physical mechanisms proposed to explain 
attenuation caused by turbulence, there have been a variety of mathematical 
models and simplifying assumptions used in order to carry out an analysis 
based on a particular assumed physical model. The different simplifying 
assumptions often lead to different end results (ref. 12). 

In an effort to bring some order to the chaotic state of turbulence 
effects on sound propagation, recourse has been made to experiment. Summar- 
ies and critiques of important experimental results can be found in refs. 

9 and 13 to 17. In addition, the books by Tatarskii (ref. 13) and Ishimaru 
(refs. 14 and 13) contain derivations and explanations of principal theore- 
tical methods used for analyzing problems of sound- and electromagnetic- 
wave propagation through a turbulent atmosphere. 

It is unfortunate, as pointed out by Brown and Clifford (ref. 9), that 
the experimental results have been complicated by extraneous factors, and 
that crucial parameters needed to validate a particular theoretical model 
were often not recorded or measured. Nevertheless, the relatively few 
experiments that have been conducted to study atmospheric turbulence effects 
have been useful in providing some indication about the relative merits of 
the competing analyses that have been conducted so far. 

It is clear that what is needed before further progress can be made 
in advancing knowledge of turbulence effects on sound propagation is a new 
experimental program designed specifically to test the premises that underlie 
the more promising physical and mathematical models. The objective of the 
study reported here is to design an experimental program — based on a spe- 
cific physical and mathematical model — which will provide all of the 
necessary Information needed to validate (or discredit) the model. The 
experiment must be practical and use state-of-the-art, obtainable equipment. 
The purpose of the experiment would be to obtain information ultimately 
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useful for determining turbulence effects on aircraft flyover noise and 
f thereby to supplement other more-fundamental research efforts. 

In the remainder of this report, we describe the specific analytical 
model chosen to represent the effects of turbulence on the att-'iuuation of 
1 sound (the model of Brown and Clifford, ref. 9), point out and describe the 

important parameters needed to employ the model and show what effect varia- 
tions in the parameters would be expected to have on attenuation of noise 
(especially from aircraft) . We then describe one analytical model for predic- 
l ting the statistical variance of the sound pressure fluctuations caused by 

atmopsheric turbulence. And finally, we describe a practical experimental 
program which would provide the necessary information to validate the 
analytical models. 

ATTENUATION OF SOUND BY TURBULENCE 
Choice of Model 

t 

To design an experiment capable of assessing the importance of atmos- 
pheric turbulence on sound propagation and providing better methods to 
predict turbulence attenuation, a working model of the physical effects 
of turbulence on sound propagation must be selected and developed, as 
required, to make the model applicable to the problem of aircraft flyover 
noise. 


The basic model chosen to predict the effect of atmospheric turbulence on 
attenuation is that developed by Brown and Clifford. In their paper (ref. 9), 
Brown and Clifford point out that the simple scattering theories such as those 
of DeLoach (ref. 15), Lighthill (ref. 19), and Blokhintzev (ref. 20) are not 
correct in that they invoke a single scatter model that does not conserve 
energy. The single-scatter models assume that the scattered acoustic energy 
is lost; they do not take into account the significant amount of energy 
scattered toward the observer (or microphone) by off-axis turbulence. 

The theory proposed by Brown and Clifford, on the other hand, is based 
on previously developed theories of forward propagation of optical waves 
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through turbulence, and appears to rest on a somewhat more solid theoretical 
foundation. The application to acoustics of results derived on the basis of 
approximations valid for optics is somewhat tenuous. However, Brown and 
Clifford argue in their paper (ref. 9} that the results they derive should 
be valid for acoustic wavelengths of up to a few meters, i.e., to relatively 
low audible- frequency sounds. 

Comparisons between calculated results based on their theory and various 
available experimental results were made by Brown and Clifford and showed 
good agreement. However, as we show later. Brown and Clifford made some 
numerical errors in their calculations. The agreement between theory and 
experiment now appears worse than reported in ref. 9. Furthermore, good 
agreement between theory and experiment can also be obtained from other 
models, such as that of DeLoach, if certain empirical parameters are chosen 
properly. 

However, the fundamental problem, mentioned in the intioductlon still 
remains — namely, that the experimental reports invariably lack key infor- 
mation on values of important parameters because the experimenters did not, 
at the time, know what parameters were most important for determining 
attenuation by turbulence. That is why new experiments based on a physical 
model are needed. 


The Brown and Clifford Model 

In this section we shall outline the derivation of the equation obtained 
by Brown and Clifford for the attenuation of sound by turbulence. In 
their model, the principal mechanism for the apparent attenuation caused 
by atmospheric turbulence is the broadening of a beam of sound. Figure 1 
illustrates a beam of sound propagating through a non- turbulent and a 
turbulent atmosphere. For a non-turbulent atmosphere, diffraction will 
cause the beam to spread from an initial diameter Dq to a larger diameter 
Df after propagating a distance L. The magnitude of depends on Dq, L, 
and the frequency f of the sound. 
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Figure 1. -Geometry of an acoustic beam propagating through a turbulent and a non-turbulent 
atmosphere. 


For a turbulent atmosphere, there is an additional spreading of the 
beam caused mainly by turbulence- induced phase fluctuations. Thus, after 
propagating a distance L through a turbulent atmosphere, the diameter of 
the beam will be equal to D which is larger than D f . The magnitude of the 
additional spreading of the beam due to turbulence depends on L, f, and 
characteristics of the turbulence. 


Neglecting atmospheric absorption (which can be accounted for separately) , 
the acoustic power carried by the sound wave is constant through any cross 
section of the beam. Therefore, the average acoustic intensity over a 
cross section is inversely proportional to the cross-sectional area and 
the following relation holds between the intensity 1^ in the absence of 
turbulence and the long-term time-averaged intensity <I> in the presence of 
turbulence: 

<I>/I f - D 2 /D 2 - D 2 /(D 2 + D 2 ) (1) 

where 7 tD 2 /4 is the additional beam area caused by turbulence- induced beam 
spreading. 

The attenuation in decibels, caused by turbulence is defined by 

A, - 10 log (I,/<I>) - 10 log [1 + (D 2 /D 2 )] (2) 

t 10 c 10 t f 
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For the quantity D 2 , Brown and Clifford use a solution obtained by 
Yura (ref. 21): 

D 2 * 25 k z/s £ J (L - s) 5/3 C 2 (s) dsj (3) 

where k * 2*rrf/c is the wavenumber of the sound, f the frequency of the sound, 
c the sound speed and s is distance along the sound propagation path from the 
source to the receiver. 


The quantity C 2 is the structure function for refractiviny fluctuations 

with units of length to the -2/3 power; it is the only parameter in the model 

that depends on the characteristics of the atmospheric turbulence. It 

will be defined and discussed in a later section. The integral in eq. (3) 

represents a sort of weighted average of C 2 over the sound path. The 

n 

integral ranges from the sound source to the receiver, i.e., s * 0 is 
the location of the source and s * L is the location of the receiver. 


For the quantity D 2 , Brown and Clifford use the following expression: 

D z = D z + (16 L 2 ) / (k 2 D z ) (4) 

In the far-field the first term in eq. (4) is negligible compared with 
the second term, and the equation is an expression of "inverse-square-law" 
spreading. The factor 16 multiplying the second term is suspect. We 
believe that it is in error by a factor of 4, i.e., the correct value is 64. 
In Appendix A we justify the value of 64. 

By combining eq. (2), (3), and (4) Brown and Clifford obtained the 
following equation for attenuation caused by turbulence: ^ 

A t - 10 log ^ ll + 1.56 k 12/s Dj^ J [(L - s)/L] s/3 C 2 (s) dsj J (5) 

L 0 

The factor 1.56 in eq. (5) is the value obtained by Brown and Clifford. If 
the factor 16 in eq. (4) were changed to 64 then the number 1.56 would be 
changed to 0.391. 



Parameters in the Brown and Clifford Model 
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The parameters appearing in the Brown and Clifford model are the 
structure function for refractivity fluctuations C*, wavenumber k, in* 
source diameter D^, and propagation distance L. In this /ir,n f we discw 
each of the parameters and illustrate what effect varyin each parameter 
would have on attenuation caused by turbulence. 

Structure Parameters * -As mentioned previously* the function C* is the 
one variable in the Brown and Clifford model that depends on the charac- 
teristics of the turbulence in the atmosphere. It can be related to two 
measurable parameters which will be defined now. The first of those Is 
the structure parameter for velocity fluctuations C*; the second is the 
structure parameter for temperature fluctuations C*. 


Before the structure parameters can be defined. It is first necessary 
to define so-called structure functions for turbulent velocities and temper- 
atures. For velocities, the structure function Is a tensor defined by the 
following expression: 


'ik 


<lv 2i - v il )(v 2k ‘ V’ 


( 6 ) 


where the symbol < > denotes a time average. The definition expresses a 
correlation of velocities, i.e., the relation of velocities at two 

-V 

neighboring points 1 and 2. The velocities at the two points are and v^. 
The subscripts i and k denote components of the velocity vectors. For 
locally isotropic turbulence depends only on the distance r between 
the two neighboring points. The component D is then just the mean- 
square relative velocity of two neighboring fluid particles along the line 
joining them. Similarly, the component is the mean-square transverse 
velocity of one particle relative to another. If the distance r between 
the two points is large compared with the size of the smallest eddies In 
the turbulent flow but small compared to the largest eddies (i.e., the 
so-called "inertial subrange" region) then it can be shown (ref. 22) that 

D tt “ (4/3) D rr and 


■ "constant" 


r 2/3 


c 2 r 2/3 

V 


(7) 
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The structure parameter for velocity fluctuations C 2 is thus defined as 
the "constant" in the "2/ 3-power law" for the structure function D^. The 
quotation marks around the word constant indicate that is actually only 
a constant over a limited region of the atmosphere and does vary, as shown 
below, with height above ground level and the condition of the atmosphere. 

Similarly, a structure function can be defined as the mean-squat e 
relative temperature difference of two neighboring fluid particles: 

D t - <(T 2 - ) 2 > (8) 

Again, it can be shown that 

D t - "constant" • r 2 / 3 » C 2 r 2 / 3 (9) 

Equation 9 defines C 2 which, like C 2 , is only a "constant" over ; limited 
region of the atmosphere. 

The structure parameters can be measured directly by, say, employing 
two sets of hot-wire anemometers located a distance r apart. However, 
for practical applications it is often necessary to relate the structure 
parameters to other quantities such as turbulent energy dissipation, which 
may be easier to measure, and then deduce numerical values of the structure 
parameters as a function of height. Practical methods for measuring 
and C 2 are discussed in a later section. 

Once C 2 and C 2 have been determined as a function of height, C 2 can be 
a v n 

found from the following equation* (ref. 13): 

C 2 - C 2 /4T 2 + C 2 /c 2 (10) 

n i v 

where T Is the ambient temperature and c is the ambient speed of sound. 

*A more exact equation for C 2 is the following: 

C 2 - C 2 /4T 2 + C 2 /c 2 + 2(0.307)C eT /4pT + (0. 307) 2 C 2 /4p 2 

where C e j is a parameter for cross correlation of humidity and temperature 
fluctuations, p is the ambient pressure, and C| is a structure parameter 
for vapor pressure fluctuations. The expanded equation was derived by 
Wesley (ref. 23) is also discussed in ref. 17. The additional parameters 
C t x and c| have not been thoroughly investigated but their effects may 
be important for sound propagation in humid atmospheres. 
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Careful measurements of the structure parameters C 2 and C 2 were made 
during two extensive experimental piogiams conducted to measure the 
characteristics of turbulence in the atmospheric boundary layer (refs. 

24 and 25). A typical set of profiles of potential temperature, wind 
speed, and wind direction, taken from ref. 24 is shown in fig. 2. The 
thickness of the convective atmospheric bo>’ndary layer is defined by 
the height of the lowest inversion base. The dashed portion of the curves 
was ootained from rawinsonde measurements which tend to be less precise 
than the balloon-borne measurements used to obtain the lower portion of 
the curves. 

By properly nondimens ionalizing the structure constants, Kaimal et at. 

in ref. 24 showed that data from all of the convective-atmosphere runs 

can be collapsed to single curves over a large range of heights. The 

nondimeasionalized data for the structure constants C 2 and cl fiom ref. 24 

v r 

are shown in fig. 3. The scaling velocity w A and temperature T* are 
defined by the following equations: 

w . * 1% Zj (g/T)] l/l (lla) 

T, ’ V”* (Ub) 

where is the surface kinematic heat flux defined by * <w’T r > at a 
height Z * 4 m; w 1 is the fluctuating component of wind in the vertical 
direction; T 1 is the fluctuating temperature; and g is the acceleration 
of gravity, (Note: we have used somewhat different notation from that 

of ref. 24). 

The dashed lines in fig. 3 are predictions for the free-convection 
portion of the atmospheric boundary layer and are based on the data 
reported in ref. 25. For heights greater than 0.1 the following 
approximate relation holds 

C 2 = constant (12) 

vi * 

Tht constant in eq. (12) is between about 1.0 and 1.5. 
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(°K) (m/sec) (degrees) 


Figure 2.*Profifes of potential temperature, wind speed, and wind direction, for 
Run 2A1. (From Kaimsi et al, ref. 24. Zj “ thickness of convective 
atmospheric boundary layer. 



Figure 3.-Vertical profiles of structure functions. (From Kaimal et al. ref. 24). 
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For the range of heights 0.01 Z <Z _< 0.1 Z^, Kalmal et al. (ref. 24) 
suggest the following approximation for the normalized C^: 

C 2 Z? /j /w 2 - 1.3 + 0.1 (Z/Z,)“ 2/s (13) 

v l w 1 

For the normalized C 2 they suggest the following approximation which Is 
Indicated by the dashed line in fig. 3: 

C 2 ZJ /j /T 2 - 2.67 (Z/Z 1 )~' ,/ 3 (14) 

Figure 3 shows that the data for the normalized C* have little scatter and 
that eq. (14) holds for heights up to about Z • 0.J 

For the typical run (Run 2A1) for which the wind and potential temper- 
ature profiles are shown in fig. 2, the following numerical values for the 
boundary layer parameters w^, T^ t and Z ^ are given in ref. 24: - 2.00 m/s, 

T a • 0.098 K, and - 1250 u. Substituting those values into eqs. (13) 
and (14) gives the following approximate expressions for C* and C^, for a 
typical midday unstable boundary layer: 

C 2 - 0.045 + 0.40 Z~ 2/ \ m'/Vs 2 (15) 

C 2 - 3.0 Z"^*, K 2 /m 2/j (16) 

where Z is in meters. 

For illustrative purposes It is useful to obtain a typical variation 
of C 2 with height. To do so, we shall assume an atmosphere with a con- 
stant temperature of T ■ 283 K. Then substituting eqs. (15) and (16) into 
eq. (10) vs obtain the following expression for a typical variation of C 2 
with height Z: 

C* - 4.0 x io~ 7 + 3.5 x 10" 6 Z -2/3 + 9.4 x ]0' s Z -,,/ \ m" 2/> (17) 

n 

The values of the coefficients in eq. (17) have been rounded to two 
significant figures. 

Brown and Clifford (ref. 9) also use the data of ref. 23 to obtain 
a typical variation of C* with height. Their equations for C* and 
are the same as eqs. (15) and (16) but the values of the coefficients 
are slightly different. If we use their values and a temperature of 283 K, 
we obtain the following equations for the structure parameters: 
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( 1 8a) 


C'/c r - 3.5 x 10- 7 + 2.9 x Ur 6 Z~ J / 3 . nT ;/3 

C‘/4T : - 9.1 x 10" 6 Z* 4/, f m~ l/3 (18b) 

and the following for the structure function C 2 

n 

C 2 - 3.5 * 10- 7 + 2.9 V ICT 6 7.~ 2 / 3 + 9.1 x 10” f Z~ 4/3 , m _;/1 (19) 
n 

In order to be consistent with Brown and Clifford we shall use eq. (19) 
instead of eq, (17) for sample calculations later. The variation of (T 
with Z* according to eq, (19)* is shown in fig. 4. Since the curve is just 
to be used for sample calculations , the range of Z has been extended 
downwards below the height for which eq. (19) is known to be A good approxi- 
mation to the variation of C* with height above ground level. 



Figure 4. -Typical variation of structure parameters with height above ground. 
Calculated according to eqs. 18 and 19. T - 283K. 
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Curves of C 2 /c 2 and C 2 /4T 2 are also shown in fig. 4. The two curves 

* cross at Z* 3.75 m. For heights greater than that C^,/4T 2 rapidly becomes 

much less than C 2 /c 2 . Near the ground, fluctuations in temperature (C^) 

dominate the value of C 2 . At heights above about 10 m, fluctuations in 

n 

wind velocity (C 2 ) control the value of C 2 . 

v n 

i 

Parameters k - V and L.-The wave number k and the sound propagation 
distance L appearing in eq. (5) are well defined and can be determined in 
a straightforward way* The remaining parameter Dq — the initial beam 
I diameter — is not well defined except for an experimental situation that 

closely approximates the beam model on which the analysis of Brown and 
Clifford is based* Thus, if an experiment were conducted using a tower- 
mounted- loudspeaker sound source, Dq would be set equal to the diameter 
^ equivalent to the area at the exit plane of the source* 

For such an experiment, the frequency of the loudspeaker-generated 
sound would be controlled through the electrical signal driving the 

loudspeaker. Thus, the wavenumber k would be Independent of diameter D n . 

» 0 

For other types of sound, such as the jet noise of an aircraft engine, 
the wavenumber k and the parameter Dq will not be independent of each other* 

The noise generation region in the jet extends aft from the tailpipe exit. 

The high frequency part of the jet-noise spectrum is created by sources 
close to the tailpipe exit and the low frequency part of the spectrum is 
caused by sources further downstream. Furthermore, the diameter of the 
jet Increases in the downstream direction. Thus the high-frequency (small 

* wavenumber) part of the spectrum is associated with a smaller characteristic 
dimension for Dq than is the low frequency (large wavenumber) part of the 
spectrum. The coupling of k and Dq can be put on a more quantitative basis 
by rewriting the factor k 12 ^ 5 D 2 in eq. (5) as (D 2 k 2 )k 2 ^ 5 . Then by intro- 

1 ducing a Strouhal number St ■ fD Q /U " W)q/ 2^M, where H » U/c is a character- 

istic Mach number, the factor k l2 / 5 D 2 can be seen to be equal to 4ir 2 M 2 St 2 k 2 ^ 5 . 
Equation (5) can therefore be written in the following form:* 

^ A c - 10 log i(j jl + 61.6 M 2 St 2 k 2/5 j^y L (a - s)/L] s /* C 2 (s) ds 

♦Again, it must be pointed out that the factor 61.6 in eq. (20) shou: 
divided by four, i.e., it should be 15.4 if the factor 16 in eq. (4) were 
changed to 64. 

» 
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Equation (20) is appropriate for application to the propagation of 
aircraft flyover noise whereas eq. (5) is applicable to noise produced by a 
loudspeaker. A comparison of the two equations shows that the frequency 
dependence of turbulence- induced sound attenuation is stronger for sound 
generated by a loudspeaker than for sound from an aircraft engine because 
of the k* 4 *^ 5 factor in eq. (5) compared with the k 2 ^ 5 factor in eq. (20). 

Furthermore, a comparison of the two equations indicates that the 
magnitude of the turbulence- induced attenuation is much less for an 
airplane-engine sound source than for a loudspeaker source. That can be 
seen by taking the ratio of the second term in the braces in eq. (20) to 
the second term in the braces in eq. (5). For small values of the second 
term the above ratio is approximately equal to the ratio of the attenuation 
for an airplane-engine sound source and a loudspeaker sound source. Calling 
the ratio R we obtain 

R - 4lt 2 M 2 St 2 /D 2 k 2 . (21) 

For reasonable values of M, Sc, and Dq (say M * Sc a 1 and Dq ■ 1 m) Che 
racio R will be much less chan unity excepc for small values of k, i.e., low 
frequencies. 

Calculation of Sound Attenuation Induced by Turbulence 
(Brown and Clifford Model) 

In this section we show the results of some calculations of turbulence- 
induced sound attenuation based on the Brown and Clifford model. 

For sound generated by a loudspeaker-type sound source, eq. (5) is 
rewritten in the following form: 

A • 10 log fl + C. k ,2 / s D 2 (L C 2 ) 6 ' 5 ) (22) 

— - ^ io [, a u nj 

where C 2 is a weighced average of C 2 over Che sound pach; i.e. 

F = (1/L) f L ((L - s) /L] s/ 3 C 2 (s) ds. (23) 

n •'O n 

The nondimens ional constant in eq. (22) is equal to 1.56 according to 
the original Brown and Clifford analysis or is eqjal to (1.56/4) ■ 0.391 
if their expression for the far-field spreading of a beam of sound in a 
non-turbulent atmosphere were incorrect by a factor of four (see Appendix A). 
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For aircraft noise eq. (20) is rewritten in the following form: 

- 10 log i 1 + C,M 2 St 2 k 2/s (LC 3 )^*} (24) 

t 10 (2 n J 

where C 2 is equal to 61.6 if the far-field beam spreading equation used 

by Brown and Clifford is correct, or to (61.6/4) ■ 15.4 if that equation 

is off by a factor of four. 

According to eq. (23), the turbulence nearer the source is weighted 

more heavily than the turbulence nearer the receiver in its effect on 

C 7 . The weighting comes about because of the [(L - s)/L] 5 ^ 3 factor in 
n 

the Integrand since s • 0 at the source and s • L at the receiver. 

2 — r 

Therefore, because C generally decreases with height (see fig. 4), C 

n « 

vould be smaller for the case of vertical sound propagation from a source 
located a distance L above a receiver than for a receiver located a distance 
L above a source. 

Figure 5 shows some geometrical relationships and defines notation 

for two cases of vertical sound propagation. One case is that of a 

sound source above a receiver, i.e., the case of an elevated source. 

The other case is that of a receiver above a sound source, i.e., the 

case of an elevated receiver. The weighted average of C 2 is given by the 

n 

following equations for the two cases: 

For an elevated source 

— * Z L / 

1C 2 - J [(Z - Zq)/(Z l - Zq)] 5 '* C 2 (Z) dZ (25) 

zo 

For an elevated receiver 

LCj -J Zl l( h - Z)/(Z L - Zo)]*/’ C 2 (Z) dZ (26) 

Z0 

Figure 6 shows calculated curves of average turbulence-induced 
attenuation coefficient in dB/m for the two cases of elevated source and 
elevated receiver. The assumed variation of turbulence structure parameter 
for refractlvity fluctuations C 2 used for the calculations is that given 
by eq. (19) and shown in fig. 4. The air temperature was assumed to be 
283 K. 
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(b) Elevated receiver. 


Figure 5.-Vertical sound propagation for a source above a receiver 
and for a receiver above a source. 



Figure 6. -Average turbulence attenuation coefficient for a vertical path. 
T = 283K, f - 4000 Hz. C, “ 1.56, a* - (1000/L) A. (dB/km). 
A t from eq. (22). 1 x 1 
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The variable In fig. 6 is the height above ground Zq of the source 
when the receiver is elevated or of the receiver when the source is 
elevated. The height Zl was always kept at 600 tn and the path length 
decreased as Zq Increased. Equations (22) and (23) with * 1.56 were 
used for the calculations. 

The calculations for fig. 6 were done for a temperature of 283 K, 
a frequency of 4000 Hz (a wavenumber k of 74.5 m” 1 ), and an initial beam 
diameter Dq * 1 m. Those values were chosen because they were the same 
as those picked by Brown and Clifford for some of their examples and they 
thus provided a check on the calculations. The quantity LC^ (i.e. ' 
the weighted average of the structure parameter C* over the propagation path 
multiplied by the path length) was found by numerical integration of 
eqs. (25) and (26). 

The results in fig. 6 indicate that when the source is near the ground 
plane and the receiver is elevated, the attenuation coefficient decreases 
by a factor of about 2 as the source is raised from 0.1 to 50 m or out 
of the region of strong turbulence. However, when the source is elevated 
at a height of 600 m and the receiver is located near the ground plane, 
then raising of the receiver does not influence the calculated attenuation. 

The difference between the calculated a c values for an elevated 
receiver and an elevated source decreased rapidly as the distance above 
the ground plane was increased, coinciding with the rapid decrease in the 

A 

C n function with increasing height in fig. 4. When there was little 
variation of C* with height above the ground plane, say for heights Z > 10 m, 
then there was little difference in the effect of turbulence on sound 
propagation. Note in fig. 6 that the values calculated for a t at 
Zq ■ 1 m are 3.8 and 5.6 dB/km. Brown and Clifford (ref. 9) calculated 
corresponding values of 1.5 and 8.0 dB/km for the same case. Their 
results were found to be incorrect due to errors in their calculations. 

Thus, the calculated difference in attenuation between the cases of an 
elevated source and an elevated receiver is much less than stated in ref. 9 
even if the value of the coefficient were equal to 1.56. 



Figure 7 shows calculated average turbulence- induced-attenuation 
coefficient for the same conditions used for the calculations displayed 
in fig. 6. However, the coefficient in eq. (22) was set equal to 0.391 
instead of 1.56. It is seen that the differences between the cases of an 
elevated source and an elevated receiver are even smaller than those shown 
in fig. 6 or reported in ref, 9. 



Figure 7.-Average turbulence attenuation coefficient for a vertical path. 

T - 283 K, f - 4000 Hz, C 1 ■ 0.391, a t ■ (1000/L)A t (dB/km) A t from eq. (22). 

Because the structure constants vary with height above the ground, the 
average turbulence attenuation coefficient for a sound path a t , in dB/km, 
varies as the sound source height is varied. Figure 8 illustrates the 
variation in a t with source height for a fixed receiver height of 1 m. 
Separate curves are shown for the two values of the constant C^. The 
values of a t shown in fig. 8 for a source height of 600 m. correspond to 
those shown in figs. 6 and 7 for an elevated source with the receiver at a 
height of 1 m. 
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Figure 8.-Average turbulence attenuation coefficient versus source height. Receiver 
is 1-m above ground level. a- t = (1000/L)A t (dB/km). T = 283 K, 

f - 4000 Hz, Dq - 1 m. A t from eq. (22). 


The variation of a t with frequency is shown in fig. 9. Separate 
curves are given for an elevated source and for an elevated receiver and for 
the two values of the coefficient C^. Note that a t is plotted on a loga- 
rithmic scale in fig. 9 in contrast to the linear scale in figs. 6 to 8. All 
four curves have the same general shape. For low frequencies the slope of 
the curves, on a logarithmic plot, is about 2.3. Thus for low frequencies 
varies as the 2.3 power of frequency. As the frequency increases the 
slope of the curves decreases. 

Figures 6 to 9 are all for a loudspeaker-type sound source, i.e., a 
sound source whose effective diameter is independent of frequency. For 
a jet-prcpelled airplane, there is a coupling between effective source 
diameter and frequency, the source size decreasing as frequency increases. 
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(dB/km) 



Figure 9.- Variation of turbulence-induced attenuation, a t , 
with frequency, f. A t according to eq. (22). 


As described previously, the coupling is through a Strouhal number relation- 
ship St ■ fOg/U. Figure 10 shows the calculated variation of o t with 
frequency, according to eq. (24) for an aircraft-type sound source. The 
characteristic Mach number M and Strouhal number St were both arbitrarily 
set equal to unity. 
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Figure lO.-Variation of turbulence-induced attenuation with frequency, f, 

for an aircraft-type sound source. M - St * 1.0. * (1000/UA^ (dB/km.) 

Aj according to eq. (24). 
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Two characteristics of the calculated curve are apparent. First, the 
variation is very closely approximated by a straight line. The slope of 
the straight line is 2/5 which implies that a t varies as the 2/5-power 
of frequency. The second important characteristic of the plot is that the 
calculated values of a t are significantly smaller than those in figs. 6 
to 9 for a loudspeaker- type sound source, thus Indicating that atmospheric 
turbulence may have little effect on aircraft noise attenuation. It 
should be kept in mind, however, that since a t is approxlmatley proportional 
to the product M 2 St? its value could Increase greatly if M or St were 
Increased from their assumed values of unity. Nevertheless, a characteris- 
tic value of unity for M and for St is reasonable and even if St were 
increased to, say, St ■ 3 the attenuation coefficient a t would increase 
by a factor of about ten and would still be rather small. 
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Calculation of Sound Attenuation Induced by Turbulence 

(Deloach Model) 

DeLoach' s equation (ref. 15) for tuxbulence-caused sound attenuation 

is based on an empirical modification of a single-scatter model. The 

equation for the attenuation over a sound path (in decibels) is 

,Z L 

A t - 1.980 k l/3 J [ (n/kLn) + sin (6 /2) ]“ s / 3 • 

Z 0 

[(C 2 /c 2 ) +0.136 (C 2 /T 2 )] dZ (27) 

where Lq is the "outer scale'* of the atmospheric turbulence and 0 c is an 
empirical parameter. 

A straightforward derivation of an equation for A t based on a single- 
scatter model would yield eq. (27) but with the equivalent angle 0^ equal 
to zero. For that case A t would be proportional to the square of the 
wavenumber. DeLoach's empirical modification leads to the introduction of 
the new parameter 6^. If sin (0 c /2) were much greater than n/kL ^ then A fc 
would be proportional to the one- third power of the wavenumber. Thus, 
DeLoach’s modification allows for a frequency dependence of A t rarging 
from a frequency-squared to a frequency to the one- third power. The exact 
frequency dependence depends on the relative magnitude of the terms r/kLQ 
and sin (0^/2) . DeLoach argues that measurements of "excess" attenuation 
of aircraft flyover noise show a one-third-power-of-the-frequency dependence. 
Also, he is able to fit experimental results for "excess" attenuation by 
choosing specific numerical values for the parameters 0 and L^. 

A deficiency of DeLoach’s model is that the quantity 0 does not seem 
to be related to any measurable physical quantity ^except for A^ itself) . 
Because 0^ cannot be measured, eq. (27) cannot be used to predict the 
magnitude of A^ even if the structure parameters C* and C* hn^e been 
measured. 

In his report, DeLoach gives several examples i ' owing that eq. (27) 
is consistent with reported measurements of excess attenuation. He defines 
excess attenuation as the attenuation a sound wave suffers that is in 
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addition to that caufed by "classical plus molecular absorption." DeLoach 
f argues that turbulence scattering is the primary cause of excess attenua- 

tion. 

Since excess attenuation is defined by DeLoach to be attenuation not 
» accounted for by viscosity and heat conduction effects,* i.e. by atmospheric 

absorption, it is important in interpreting experimental data to accurately 
account for the effects of absorption. That fact can be illustrated by 
looking at some data that DeLoach used as a check of his theory. The data 
^ vere taken from ref. 26 and represent a composite of measured total atten- 

uation coefficients from several experiments. Temperature and humidity 
were measured during the experiments but the structure parameters C* and C 
were not. The outer scale of turbulence Lq was rot measured either. 

» 

Figure 11, taken from ref. 15, shows the measured data along with 

three calculated curves. Two of the curves were calculated by DeLoach, 

the third was calculated as part of this study. The short-dashed straight 

^ lines show the atmospheric absorption calculated by DeLoach according to 

a method described in ref. 15. Tba solid-line curves show the sum of the 

calculated absorption and tie attenuation caused by scattering calculated 

according to eq. (27). The numerical values shown for the structure 

parameters were assumed values. The parameters 6 and L. were determined 

> cu 

by DeLoach by requiring that eq. (27) give a best fit to the data accord- 
ing to the method of least squares. 

It can be seen from fig. 11 that the data can be closely matched by 

* a proper choice of 0 and L ft . On the other hand, the measured attenuation 

c u 

can also be explained on the basis of absorption alone provided that the 
absorption is calculated according to another method. The long-dashed/ 
short-dashed curves show absorption calculated according to the method of 

* ref. 3. The procedure of ref. 3 for calculating pure- tone atmospheric 
absorption coefficients is based on a more up-to-date understanding of the 
physical phenomena than the pro .dure used by DeLoach in ref. 15. 

^ *Viscosity effects include effects associated with both of the two coeffi- 

cients of viscosity (e.g., see ref. 27) that are in the Navier-Stokes equa- 
tions. Therefore, both "classical" absorption and molecular relaxation 
effects are included here as effects of viscosity. 
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DATA 

ABSORPTION ALONE ACCORDING TO REF. 16 

SUM OF ABSORPTION AND SCATTERING ACCORDING TO REF. 16 

ABSORPTION ALONE ACCORDING TO REF. 3 
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ABSORPTION ALONE ACCORDING TO REF. 15 

SUM OF ABSORPTION AND SCATTERING ACCORDING TO REF. 15 

ABSORPTION ALONE ACCORDING TO REF. 3 
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ABSORPTION ALONE ACCORDING TO REF. 15 

SUM OF ABSORPTION AND SCATTERING ACCORDING TO REF. 15 

ABSORPTION ALONE ACCORDING TO REF. 3 
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C y - 0.15 m 2/3 /**c; C T - 0.17° C/m 1/3 ; 
0 C » 0.67°; Lq- 37.7 m;T-12.1°C; 
and ralativa humidity * €8.3 ptroant; 



TEMPORAL FLUCTUATIONS CAUSED BY TURBULENCE 


Turbulence in the atmosphere has two effects on sound measured in the 
far field of a noise source* The apparent reduction of the mean intensity 
of the sound has been discussed in the previous section. In this section 
ve discuss the effect of atmospheric turbulence on the variation of 
sound intensity with time during a measurement. 

Large fluctuations in the amplitude of the sound pressure level have 
been known from common experience and from outdoor sound propagation experi- 
ments for some time 9 though the cause has been a matter of speculation. 

In ref. 28, Knudsen describes the results of an experiment conducted 
on the campus of UCLA in 1934. A transmitter and a receiver were set up 
about 100 ft apart. The experiment was conducted early on a Sunday morning, 
before sunrise, with no perceptible wind. With a 4000-Hz tone as a signal, 
"the level at the receiver fluctuated violently over a range of more than 
10 d& with short periods (0.1 s or less) and long periods (several seconds) 
all jumbled together." In an experiment conducted on a Mojave dry lake 
in the summer of 1940, the fluctuations were found to increase with sound 
frequency and distance between transmitter and receiver. The major cause 
of the fluctuations was found to be inhomogeneities in the temperature 
and velocity of the atmosphere. 

Rudnick (ref. 29), following up the observations of Knudsen and Delsasso, 
investigated the propagation of sound waves in an anechoic chamber wherein 
the air temperature was constant at 20°C. He introduced temperature inhomo- 
geneities by heating a 5-m-long resistance wire stretched across the room. 

He found that the fluctuations in received sound pressure level varied with 
frequency and angle of incidence of the sound wave onto the sheet of hot 
air rising from the heated wire, the largest fluctuations occurring near 
grazing incidence and at higher frequencies. A characteristic of the 
fluctuations was the presence of short period fluctuations — of the order 
of seconds and fractions of seconds — superimposed on variations with 
periods of the order of minutes. For a given angle of incidence, the 
fluctuations at different frequencies could be correlated statistically by 



the ratio of the root-mean- square deviation of the sound pressure about 
a given amplitude to the square of the frequency of the sound. 

As part of a Symposium on Aircraft Noise in November 1952, Ingard 
(ref. 30) presented a review of meteorological effects on sound propaga- 
tion including a discussion of the effects of turbulence. Velocity 
fluctuations associated with the gustiness of the wind were shown to 
produce fluctuations that increased with gustiness, frequency of the sound, 
and distance between transmitter and receiver. The fluctuations ranged 
from 10 to 20 dB at 2000 Hz and from 5 to 25 dB at 4000 Hz for wind speeds 
from 6 to 11 m/s. The source of the fluctuations was thought to be 
atmospheric turbulence. 

Wiener and Keast (ref. 31) reported the results of an extensive study 
of the effects of atmospheric conditions on sound propagation outdoors 
from loudspeaker sound sources. At night when there was a strong positive 
temperature gradient (or inversion) because of radiation cooling of the 
ground, the air near the ground was stable and turbulent wind fluctuations 
were limited in amplitude and of low frequency. Under these conditions 
the fluctuations in the sound pressure level at the microphones were 
small. On sunny, windy days, the air temperature profile showed a large 
negative gradient (lapse), the turbulence was strong, and the air near the 
ground was unstable. Large fluctuations in the received sound pressure 
levels were noted. The fluctuations contained appreciable high-frequency 
components. 

Defining peak-to-peak fluctuations as the difference between the 
maximum and the minimum sound pressure levels in a 30-s observation period 
and plotting che results as a function of distance from the sound source 
for various octave band center frequencies gave a measure of the effect 
of atmospheric turbulence on the fluctuations of the received sound pressure 
level. F^r typical unstable daytime conditions and upwind propagation, 
the peak-to-peak fluctuation at 200 ft varied from 12 dB for the 425-Hz 
band to 25 dB for the 3400-Hz band; at 3200 and 4500 ft Lhe fluctuations 
ranged from 5 to 7 dB. For typical stable nighttime conditions and propa- 
gation downwind, there was less variation with distance and the fluctuations 


28 



I 


ranged from 3 to 5 dB for the 425-Hz octave band to 5 to 13 dB for the 3400- 
t Hz octave band, with the larger values occurring at the greatest distances. 

The periods of the fluctuations were generally short for either atmospheric 
condition. 

I In a study of the propagation of sound over ocean waters in fog, 

Wiener (ref. 32) found little effect of atmospheric turbulence on sound 
propagation. This result was attributed to the very stable condition of 
the atmospheric boundary layer over water. 

» 

Wiener, Malme, and Gogos (ref. 33) describe the results of experiments 
on the propagation of sound in city streets. They found that the temporal 
fluctuations in the received sound pressure level were less than noted in 
I open country but were still significant. At night when the atmosphere 

was reasonably stable, the maximum peak-to-peak fluctuation was about 3 dB 
and increased with signal frequency and d .stance from the sound source. 
Typical day-time fluctuations were about 8 to 10 dB. The period of the 
I largest fluctuations was on the order of a few seconds. 

All of the experiments on sound fluctuations that we have just 
summarized suffer from the same basic defect that existed for the experi- 
^ ments that were conducted to measure attenuation caused by turbulence. 

That is, because of a lack of a physical model, important parameters that 
influence the experimental results were not measured. 

We shall now give some results, based on weak- fluctuation theory 
which should provide guidance for experiments to determine atmospheric 
turbulence effects on sound fluctuations from either a static sound source 
or an airplane during a flyover noise test. 

* A statistical measure of the fluctuations in intensity is the varince 

of the intensity, i.e., the mean-square of the fluctuations in the intensity 
of a sound wave relative to the mean intensity. 

’ For weak fluctuations, a simple solution exists for the variance of 

the logarithm of the Intensity (e.g., see refs. 14 and 17) 
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oJ nI = <(lnl - <lnl> ) 2 > - 4o 2 

* 4 (0.5628) k 7/6 L ll/6 [C 7 ] 

n x 

where a 2 is the "log amplitude variance" and [C 5 *] 

X n X 

wave than for a spherical wave. 

For a plane wave: 


(28) 

is different for a plane 


[F} x * (l/L) f L (1 - s/L ) 5/5 C 2 (s)ds 


(29) 


For a spherical wave: 

[C 2 ] » (l/L) r L (s/L) s/6 (1 - s/L) s/6 C 2 (s)ds (30) 

n X +q n 

For a 2 «1 it can be shown (ref. 4) that 
X 

o 2 * [<I 2 > - <I> 2 ]/d> 2 » exp a 2 nI - 1 (31) 

The above equations were originally derived on the basis of approximations 
valid for optical wavelengths. The same is true for the equations used 
as the basis of the "excess attenuation" model of Brown and Clifford. 
However, as mentioned previously, Brown and Clifford have shown that their 
results are also valid for wavelengths appropriate to sound propagation 
at audible frequencies. It is likely that the equations for the variance 
of intensity are also applicable to acoustic wave propagation problems. 


Inspection of eqs. (28*31) shows that the same parameters appear as in 

the Brown and Clifford equations for turbulence- induced sound attenuation. 

Those parameters are the wavenumber k, the propagation path length L, 

and the structure parameter for refractivity fluctuations C*. On the basis 

of the discussion in the previous section it would be anticipated that the 

variance of the intensity would therefore be largest when C 2 is largest, 

n 

i.e., near the surface of the ground. This result agrees with observa- 
tions from field measurements of a variety of noise sources ranging from 
aircraft to motor vehicles and railroad r.rains. 

EXPERIMENTAL VERIFICATION PROGRAM 
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There are several reasons why reports on previous experiments lacked 
information on key parameters: First, some experiments were carried out 

without having a theoretical model available when the tests were planned. 
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Consequently, it was not clear as to what parameters were Important to 
measure. Second, many experiments which were later used to study atmospheric 
turbulence effects were conducted primarily for some other purpose. There- 
fore, atmospheric turbulence characteristics were not measured. That was 
usually the case for flyover noise measurements, for example. And finally, 
at the time that some of the experiments were conducted, there were no 
readily available, practical procedures for measuring the necessary turbu- 
lence parameters. 

In this section of the report we describe two general experiments 
designed to validate (or discredit) the Brown and Clifford model. The 
experiments would also be useful for testing other models of sound propa- 
gation through a turbulent atmosphere since other models, such as DeLoach's, 
contain most of the same parameters that are in the Brown and Clifford 
model . 

The two general types of experiments are: (1) tower experiments and 
(2) aircraft flyover noise experiments. For a tower experiment, a loud- 
speaker is mounted on top of a tall tower. Meteorological data are obtained 
from instruments mounted to the tower at various heights. Acoustical data 
are obtained from microphones on the ground and attached to the guy wires 
that stabilize the tower. Such a tower experiment is described in ref. 34. 
The tower used was a 150-m tower near Haswell, Colorado. Although part of 
the experimental program involved a study of sound pressure fluctuations 
caused by turbulence, structure parameters were not measured. 

Numerous other large towers exist and have been used extensively to 
study the atmospheric boundary layer. In fact, in his review article on 
the atmospheric boundary layer, Panofsky (ref. 35) calls the lower part 
of the boundary layer the tower layer because of the numerous tower observa- 
tions that have been conducted to study that region of the atmosphere. 

An example of a tall tower that would be suitable for a study of 
atmospheric turbulence effects on sound propagation is the 300-m instru- 
mented research tower near Boulder. Colorado. The tower is operated by the 
Wave Propagation Laboratory of the National Oceanographic and Atmospheric 
Administration (NOAA) . 
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Aircraft flyover noise experiments can be conducted with an aircraft 
climbing, descending, or flying level. However, level-flight flyovers 
are easiest to interpret. The aircraft engine (or engines) itself is the 
primary sound source. Microphones are located near the ground beneath 
the flight path of the airplane. For a study of atmospheric turbulence 
effects the easiest data to analyze and interpret would be those data 
obtained from microphones directly beneath the flightpath and taken at 
times when the airplane is nearly overhead. For those conditions refraction 
by steady wind or temperature gradients and ground attenuation effects 
are minimized. 

Meteorological data for aircraft flyover noise tests are ordinarily 
obtained from Instruments carried aboard a light airplane. The meteoro- 
logical or "met" airplane generally obtains data during a spiraling descent 
from a height above the greatest height of the test airplane to a height 
of 5 to 10 m above the ground. Meteorological data from the met airplane 
should be supplemented by data from a ground-based weather station located 
in the vicinity of the microphones. 

An experiment to verify a model for turbulence-induced sound attenua- 
tion should measure the important parameters that enter into the model. 

For the Brown and Clifford model they are the wavenumber k of the sound, 
the sound-source diameter Dq, the propagation distance L, the structure 
parameters C* and C*, and the temperature T. The sound speed c can be 
deduced from the temperature; and the wavenumber spectrum can be deduced 
from the frequency spectrum and speed of sound. 

The basic procedure for experimental verification of the model would 
be to measure the sound pressure level spectrum at microphones located at 
fixed locations for a range of conditions of atmospheric turbulence and 
sound propagation path lengths. The range of atmospheric conditions 
would mean a range of the structure parameters C* and C* as well as a 
range of air temperature and humidity. Corresponding to the different 
atmospheric conditions would be different measured and predicted attenua- 
tions caused by atmospheric turbulence. 
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It is also important to measure other atmospheric properties in addition 
to those that influence attenuation caused by turbulence. That is required 
because it is unlikely that, from test to test, only the turbulence 
characteristics of the atmosphere would change. Other quantities that 
determine the amount of atmospheric absorption of sound will also vary. 

Since absorption is undoubtedly the principal sound attenuation mechanism 
(in addition to geometric spreading, of course) it is important to accurately 
adjust the data to be able to distinguish between the attenuation caused 
by absorption and that caused by turbulence effects. The primary atmospher- 
ic characteristics that determine the amount of absorption are temperature 
and humidity. Pressure has a relatively small effect. Therefore, it is 
necessary that accurate measurements be made of temperature and humidity 
over the entire lengths of the sound propagation paths. A measurement 
of barometric pressure should be obtained at the surface station near 
the microphones. 

Since temperature and humidity measurements are routinely made and 
procedures for making the measurements are well established, we need not 
delve further into methods for their measurement here. Measurements of 
the structure parameters C 2 and C 2 , however, are not routinely made for 
sound-propagation experiments, and since they are important parameters 
influencing sound attenuation by turbulence we shall now describe methods 
for their measurement. 

Measurement of C* 

Although it is possible to measure directly by correlating the 
signals from two hot-wire anemometers located a fixed distance apart and 
making use of eq. (7), it is much easier to obtain C 2 by an Indirect 
method. The indirect method depends on the following relation between C 2 
and the dissipation rate of the kinetic energy c: 

C 2 - Aaie 2 / 3 (32) 

A derivationof eq. (32) can be found in ref. 36. The quantity oti is a 
nondimensional universal constant which has been found to be about 0.5 
from a variety of independent measurements. 
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There are several practical ways to measure e and thus, from eq. (32), 
to determine C 2 . A rather direct measurement of z has been developed by 
Wyngaard and Cot£ (ref. 37). They measure the fluctuating velocity u(t) with 
a single hot-wire anemometer. Next, u(t) is differentiated to obtain 3u/3t. 
Taylor’s frozen turbulence hypothesis is then invoked to obtain 3u/3x. 
Finally, e is obtained from the relation e * 15v <( 3 u/ 3 x) 2 > where v is the 
kinematic viscosity of the air (see ref. 22 , pp. 123-126). 

The "direct" measurement of e is not the simplest way to determine the 
dissipation. Its primary use has been, rather, for determining the value 
of the constant 04 that appears in eq. (32). A simpler way to determine z 
is from measurements of the frequency spectrum of the u-component of 
atmospheric turbulence. The u-spectrum of turbulence can be determined by 
measuring the velocity component u(t) with a sonic anemometer and using a 
fast-Fourier-transform technique to compute the spectrum. In the inertial- 
subrange region, Kolmogorov’s law for the u-spectrum is 

F u (k l ) = a l e2/3 k I 5/3 (33) 

where represents the wavenumber spectrum of the turbulence and 04 Is the 
same constant that appears in eq. (32). 

Taylor’s hypothesis of frozen turbulence can be used to convert eq. (33) 
from an expression for the u-spectrum in wavenumber space to an expression 
in frequency space. According to Taylor’s hypothesis the following relation 
holds between wavenumber k^, frequency f^, and mean wind speed U^: 

k L - 2TTf 1 /U Q (34) 

If we denote the u-spectrum in frequency space by the symbol ( f ^ , 
the following equation holds for the mean square of u(t): 

/ CD ~ 00 

F u (k 1 )dk 1 - <u 2 > - J S u (f 1 )df 1 ( 35 ) 

0 0 

Therefore, with dk^ ■ (k^/f^)df^ from eq. (34) we require that 

HW ■ £ i W < 36 > 

Substituting eqs. (34) and (36) into eq. (33), the following equation is 
obtained for the u-spectrum in frequency space: 
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S u (f 1 ) - (2ir/U Q )- 2 / 3 a ie 2 / 3 f " s/3 (37) 

Equation (37) can be solved for e In the following form: 

e l/3 ■ (2ir/0 Q ) 1 / 3 «i l/2 fj /s [S u (f 1 )] l/2 (38) 

Equation (38) Is the final expression needed to determine e from a 
measured spectrum of S^(f^) vs fj^ For example, for f^ - 5 Hz there would 
be a corresponding measured value of S u (f^) of S u (f^) • S u (5). Then eq. (38) 
would yield a value of e 1 ^ 3 of 

e 1 / 3 - 2tt/U 0 ) i/3 ai 1 ^ 2 5 s/s [S u (5 )] l/2 (39) 

which can be solved knowing the value of Uq and 04 . 

The methods just described for obtaining e rely on Instruments 
that would be difficult to use In a "met" airplane. Therefore, they might 
be considered for a tower- type experiment. Fortunately, instruments 
have been developed, and are commercially available, for measuring e In 
an airborne system. The system is called the Universal Indicated Turbulence 
System (UITS) by its manufacturer Meteorological Research, Inc. in 
Altadena, California. 

The design of UITS is based on eq. (37). A pressure sensor senses 
the fluctuating pressure from a pitot tube on the airplane. The electrical 
output voltage from the pressure senior is proportional to pU 2 , where 
U ■ Uq + u(t) with Uq the mean true airspeed of the airplane and u(t) 
the fluctuation of the true airspeed about the mean. Therefore, the mean 
output voltage from the sensor is proportional to pU 2 and the fluctuating 
output voltage is proportional to 2pU^u + pU 2 * 2pU Q u for small values of 
u(t). 


The power spectrum of the fluctuating output voltage is thus propor- 
tional to (PUq ) 2 times the power spectrum of the fluctuating velocity u. 

If we denote the power spectrum of the sensor output voltage by C(f^), then 

G(f L ) « (pU Q ) 2 S u (f x ) (40) 

where S u (f^) is given by eq. (37). 
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The electrical power at the output of the UITS system is proportional 
to a quantity R 2 defined by 

R 2 •J G(f 1 )B 2 (f 1 )df 1 (41) 

0 

where B 2 (f^) is the frequency response of the filter through which the 
fluctuating voltage from the pressure sensor is passed. The filter is a 
low-pass filter designed so that 3 2 (f^) is approximately zero for frequen- 
cies outside the inertial subrange and approximately unity for frequencies 
within the inertial subrange. 

Using eqs, (41), (40), and (37) the following equation is obtained for 
the root-mean-square value of R 

R - ✓S 7 = constant (pU 2 ) 2 / 3 (pe/p ^) 1 / 3 (42) 

where P re f is a constant reference density. 

From eq. (42) it is seen that the output of the "epsilon-meter" should 
be proportional to (pU 2 ) 2/3 (pe/p rgf ) l / 3 . However, in the UITS equipment, 
a special electrical circuit is provided which alters the gain of the 
measuring system by a factor of (pU 2 )" 2 / 3 . Thus R is proportional to 
(P £ /P re f ) X ^ 3 independent of the airplane speed. The constant of proportion- 
ality must be obtained by calibrating the system prior to the test. Then, 
knowing R, air temperature, air pressure, and P ref > the value of e ly ^ 3 can 
be determined and hence C 2 using eq. (32) with 04 * 0.5. The ratio C 2 /c 2 
can then be found for the values of c calculated from the temperature. A 
more complete description of the UITS design can be found in ref. 38. 

Measurements of e have been made, using UITS, for some flyover noise 
tests. But little use has been made of the data for atmospheric propaga- 
tion studies. For example, UITS equipment was used for the DC-9 Refan 
flight demonstration program (ref. 39 ), but only a small sample of the 
measured e data was reported, and no attempt was made to determine values of 
the structure parameters from the measured dissipation data. The full 
set of recorded data represents a valuable store of useful information for 
a study of the effects of atmospheric turbulence on sound propagation since 
both meteorological data — including profiles of e — and sound pressure 
level data were obtained during the flyover noise tests. 


36 



VBp, . 


t 


t 


» 


I 


ft 


I 


t 


I 


» 


I 


I 

Another flight test that was made, using UI7 , is described in ref. 40. 
The purpose of the test was to evaluate the UITS system by comparing values 
of e from the UITS system against values of e measured from a different 
type of instrument that was mounted on a tower. The tests were accom- 
plished by flying the airplane carrying the UITS equipment past the tower. 
The results were that the values of e, determined by the two instruments, 
agreed with each other within + 10 percent. 

Measurement of Cy 

It is possible to determine the structure function C 2 in an indirect 
way, analogous to that used to obtain C 2 . The method depends on the 
following equation which is derived, for example, in ref. 36: 

C 2 - 4BiNe" 1 / 3 (43) 

-'where N is the dissipation rate of <T' 2 /2> and T' is the fluctuating 
temperature. The factor Si Is a universal constant whose value has been 
determined to be close to 6i * 0.8. Values of N can be determined from 
measurements of other quantities by making use of the "budget" of <T ,2 /2>: 

0 - - <w'T’> (3T/3Z) - (1/2) 0<w'T'>/3Z)-N (44) 

Below about 50 m from the surface, the second term on the right hand side 
of eq. (44) is negligible and N can be determined within about 10 percent 
from <w , T’> (3T/3Z) (ref. 35). 

Although it is useful to use such an indirect method to determine C^, 
it is usually easier to measure C 2 directly, at least for a tower-type of 
experiment. Equipment is available commercially for the direct measurement. 
One manufacturer of a sensor is Atmospheric Instrumentation Research 
Company (A.I.R. Co.) of Boulder, Colorado. The A.I.R. C^ sensor utilizes 
two rapid-response temperature probes separated a fixed distance apart. 

The structure parameter C 2 is determined from the relations given in eqs. 

(9) and (10). The A.I.R. type of C^ sensor has been used often for measure- 
ments from towers and from tethered balloons. It could not be directly 
used to measure C 2 from a "met" airplane. However, since C 2 rapidly 
becomes small compared to C 2 with increasing height above the ground 
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(as shown in fig. 4) , it would ordinarily not be necessary to measure 
vertical profiles of C* for a flyover noise experiment, though it would be 
desirable to measure C* at locations near the ground. 

CONCLUSIONS AND RECOMMENDATIONS 

1. In any experiment to evaluate the effects of atmospheric turbulence 
on sound propagation it is crucial to recognize and to measure the important 
parameters that influence the results. This study has shown that the two 
Important parameters, for sound propagation, that characterize atmospheric 
turbulence are the structure parameters C* and C*. Those parameters can 
be easured in regions of the atmosphere where they are important by using 
practical state-of-the-art equipment and procedures described in this 
report. 

-• It is also important to make measurements of other parameters that 
are needed to adjust the data for extraneous factors. For a study of the 
apparent attenuation of sound caused by turbulence it is Important to 
measure profiles of temperature and humidity over the sound propagation 
path so as to be able to accurately account for the effects of sound 
absorption by the atmosphere. It is also important to use an accurate 
method to calculate the sound absorption. The method described in ref. 3 
is recommended. 

3. The model for apparent attenuation of sound by turbulence developed 
in this study predicts that the magnitude of the attenuation will be less 
from a flyover noise experiment than for a tower experiment utilizing a 
stationary sound source. The analysis also predicts a much less rapid 
increase in attenuation with frequency for a flyover noise test than for a 
tower test. Both effects are caused by the fact that for a sound source 
such as the jet of an aircraft engine there is a coupling between the wave- 
number of the sound and the effective source diameter. In contrast, for a 
loudspeaker type of sound source the wavenumber and source diameter can be 
varied independently of each other within limits. 


38 



4. The same parameters that are important for attenuation of sound by 
turbulence, namely C* and C*, are also the important parameters for fluc- 
tuations in sound Intensity caused by turbulence. Therefore, sound fluctu- 
ations can also be studied as part of an experiment designed to study the 
apparent attenuation of sound caused by turbulence. Because of the transient 
nature of aircraft flyover noise, it is always desirable to use several 
microphones and to have several repeat runs in order to obtain an ensemble 
average which will increase the statistical confidence in the results. 

For a study of the variance of the fluctuations in the intensity of air- 
craft noise, ensemble averaging would be a necessity in order to calculate 
the statistical quantities needed. 

5. Aircraft flyover noise data now exist (ref. 39) which Include 
measurements of the important parameters that determine turbulence effects 
on aircraft sound propagation. The data have not been used to study the 
turbulence effects. The existing data should be analyzed before further 
flyover noise tests to study sound propagation are conducted. 
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APPENDIX A 

THE FAR-FIELD DIAMETER OF A BEAM WAVE 

Brown and Clifford (ref. 6) give Che following expression for che 
diameter of a beam of sound propagation in a non-turbulent atmosphere 
[tq. (4) ] 

Dj « D 2 + (16 L 2 )/(k 2 D Q ) (Al) 

The purpose of this Appendix is to show that the number 16 in the second 
term is probably too small by a factor of 4, i.e. , that the correct number 
is 64. 

First of all, eq. (A ) with the number 64 replacing the number 16 agrees 
with a result derived by Ishimaru (ref. 14) for a beam wave in free-space. 

He derived an approximate solution for the case of a beam wave having a 
Gaussian amplitude distribution and a parabolic phase distribution with 
a radius of curvature R c , at the source location x * 0 [see fig. Al]. 
Ishlmaru's equation for the beam diameter at a distance L from the source 
is a function of a complex parameter a and is 

D 2 * Dq [(1 - C^L) 2 + (O.L) 2 ] (A2) 

where a * a + ia. - (4A/xD 2 ) + i(l/R ). 

r x u c 

Substituting the expression for a into eq. A2 we obtain the following 
equation: 

D 2 - D 2 (1 - L/R c ) 2 + 16A 2 L 2 /(tt 2 D 2 ) 

- D 2 (1 - L/R c ) 2 + (64L 2 /k 2 Dj) (A3) 

For the collimated buam case, where R becomes infinite, eq. (A3) agrees 
with the expression used by Brown and Clifford in eq. (Al) except that 
Brown and Clifford's factor of 16 is replaced by the factor of 64 in 
Ishlmaru's equation. 

As a further check, consider the problem of an acoustic source 
consisting of a piston oscillating in an infinite plane baffle (fig. A2) . 
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Figure A1. -Geometry of acoustic beams having Gaussian amplitude distributions and parabolic 
phase distributions with a radius of curvature R at the source (after Ishimaru, 
ref. 14). 

That problem has been considered in ref. 41 where the following asymptotic 
solution for the velocity potential is givens 

tfi - (Va 2 /r)e ilcr " iut (J^ka sin 6)/ka sin 6] (A4) 

where the piston is oscillating with a velocity Ve” iwt , a « D Q /2 is th; 
radius of the piston, oj * 2iff, and is a Bessel function of the first 
kind of order 1. The solution is valid when the distance r is much greater 
than the radius a. 
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Figure A2.-Piston oscillating in its own plane, 
piston radius * Dq/2 . 


f From eq. (A4) the acoustic power per unit area can be determined to 

be proportional to 

q(r, 8) ** [ajpkV 2 a 4 / (2r 2 ) ] [J^(ka sin 9)/ka sin 0] A (A5) 

^ where p is the ambient air density. 

The beam diameter D is defined by the relation q(r, 0g)/q(r, 0} * 1/e, 
where D * 2r sin 0 q. From eq. (A5) 

I (Jj^ka sin 0 Q )/ka sin 6 Q ] 2 /(l/4) - 1/e (A6) 

Equation A6 is satisfied for ka sin 0^ ■ 1.915. Therefore, the beam 
diameter is given by 

D - 2r sin 0 Q - 7.660 r/kD Q (A7) 
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For narrow beams, r * L. Therefore, from eq. (A7) we obtain the 
following equation for the square of the beam diameter in the far field 

D 2 - 58.68 L 2 /k 2 D 2 (A8) 

Comparing eq. (A8) with eq. (Al) for the far field shows that the factor 
multiplying the second term should probably be closer to 64 than to 16. 
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APPENDIX 8 


SYMBOLS 

a piston radius Dq/ 2 (fig. A2) , m 

A t apparent attenuation caused by turbulence, dB 

B z (f^> relative frequency response of filter used in U1TS system 

c speed of sound, ms'" 1 

C 2 structure parameter for vapor pressure fluctuations, N 2 m -1 4 ^ 3 

C _ parameter for cross correlation of humidity and temperature 

e fluctuations, KNm -8 ' 3 

C 2 (s) structure parameter for refractivity fluctuations, m ” 2 / 3 
C 2 weighted average of C 2 over a sound path (eq. 23), m -2 / 3 

[C 2 ] weighted average of C 2 over a sound path (eqs. 29 and 30), m" 

n x ^ 

C 2 structure parameter for temperature fluctuations, I^m** 2 / 3 

C 2 structure parameter for velocity fluctuations, m 4/ ^ 3 s" 2 

nondimensional constant (eq. 22) 
nondimensional constant (eq. 24) 
beam c. ameter in absence of turbulence (fig. 1), m 
ik-component of velocity structure function tensor, m 2 s“ 2 

D fr rr-component of velocity structure function tensor, m 2 ^ -2 

D^, structure function for temper ature, K 2 

increase in beam diameter caused by turbulence, m 
D Q initial beam diameter (fig. 1), m 

F y (k^) u-spectrum of turbulence in wavenumber space, m 3 s“ 2 
f frequency of sound wave, Hz 

f^ frequency of turbulent velocity fluctuations, Hz 

g acceleration of gravity, ms” 2 

I acoustic intensity, Win -2 
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*2 

T « 


U 

u(t) 

U 0 

V 


v li* v lk 

V 2i' V 2k 
w' 


long-term term-averaged acoustic Intensity, Wm -2 

acoustic intensity in the absence of turbulence, Wm -2 

Bessel function of the first kind of order 1 

acoustic wave number, m -1 

sound propagation distance, m 

outer scale of atmospheric turbulence, m 

characteristic Mach number 

dissipation rate of <T ,2 /2>, K 2 s -1 

ambient pressure, Nm~ 2 

<w'T'> at Z ■ 4m, surface kinematic heat flux. Kins -1 
distance between two fluid particles, m 
radius of curvature of acoustic beam (fig. Al) , m 
distance along soundpath from source to receiver, m 
Strouhal number 

u-spectrum of turbulence in frequency space, m 2 s -1 

ambient temperature, K 

temperature at "point 1", K 

temperature at "point 2", K 

fluctuating temperature, K 

scaling temperature (eq. lib), K 

Uq + u(t), true airspeed, ms -1 

fluctuating component of true airspeed, ms -1 

mean true airspeed, ms -1 

amplitude of piston velocity, ms -1 

i'th and k'th velocity component, respectively, at point 1, ms -1 
i'th and k'th velocity component, respectively, at point 2, ms -1 
fluctuating component of wind in vertical direction, ms* 1 
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scaling velocity (eq . 11a) , ms” 1 
Z height above ground, m 

thickness of the convective atmospheric boundary layer, m 
Z^ height of elevated source or receiver (fig. 5) , m 

Zq height of source or receiver (fig. 5), m 

a complex parameter in equation for beam diameter, (eq. A2), m~ 1 

a 1 imaginary part of a, m T 1 

a r real part of o, m _l 

a t (1000/L)A £ , average turbulence attenuation coefficient, dB m -1 

nondimens ional constant (eq. 32) 

8^ nondimens ional constant (eq. 43) 

£ dissipation rate of kinetic energy, m 2 s” 3 

8 angular coordinate (fig. A2), rad 

0 c empirical parameter (eq. 27), rad 

0q value of 8 at "edge" of acoustic beam 

\ acoustic wavelength, m 

v kinematic viscosity, m 2 s~ x 

p ambient density, kg m~ 3 

P re j reference density, kg m~ 3 

o 2 normalized variance of acoustic Intensity (eq. 31) 

o 2 ^^ variance of the logarithm of acoustic intensity (eq. 28) 

o 2 log amplitude variance (eq. 28) 

A 

\p velocity potential, m 2 s~* 

u> 2irf, angular frequency of oscillating piston, rad s'* 1 
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